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Abstract 



Non-equilibrium excess noise of a short quasi one-dimensional constriction 
between two superconductors is considered. A general expression for the 
current-current correlation function valid for arbitrary temperatures and bias 
voltages is derived. This formalism is applied to a current-carrying quan- 
tum channel with perfect transparency. Contrary to a transparent channel 
separating two normal conductors, a weak link between two superconductors 
exhibits a finite level of noise. The source of noise is fractional Andreev scat- 
tering of quasiparticles with energies \E\ greater than the half- width A of the 
superconducting gap. For high bias voltages, V> A/e, the relation between 
the zero-frequency limit of the noise spectrum, S(0), and the excess current 
/exc reads S(0) = (l/5)|e|I ex c- As A — > both the excess noise and the excess 
current vanish linearly in A, their ratio being constant. 
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Introduction. As is well known ((see, e.g., Ref. 0), electron transport through short su- 
perconducting weak links can be described in terms of multiple Andreev reflections (MAR). 
Calculations of dc / — V curves of voltage-biased weak links have predicted specific 
current singularities at voltages V = 2 A/en, n = 1,2, where 2 A is the superconducting 
energy gap and e the electronic charge. This feature is a manifestation of MAR. However, 
dc properties carry only indirect information about weak link dynamics - more direct infor- 
mation can be found from theoretical studies of ac currents. Studies of ac transport have 
mostly been carried out for large voltages, eV 3> A, the limitation being caused by the fact 
that at low voltages it is necessary to take into account a large number of MAR. Recently, 
Averin and Bardas || have considered a model of a short constriction between two super- 
conductors which permits quantitative calculations of ac currents for arbitrary voltages V 
and transparencies D. Unfortunately, experimental studies of ac characteristics are difficult. 

In addition to studies of transport properties, investigations of non-equilibrium excess 
noise provide qualitative information about the dynamics of mesoscopic junctions. This type 
of noise can be investigated without ac bias, and directly probes the correlation between 
particles passing through the junction. It is not surprising that noise properties of weak links 
of the Josephson type have been extensively discussed (see, e.g., Refs @.|1(J for a review). 
The specific role of quantum mechanical correlations in normal contacts due to the Fermi 



statistics of electrons was considered in Refs. \\ \ -13]. There it was found a specific quantum 



reduction of noise in comparison with the so-called Schottky limit |L4 . 

In the present study we discuss noise properties of a short voltage-biased constriction 
separating two superconductors. The junction is modeled as an adiabatically-smooth point 
contact with typical length L much larger than the Fermi wave length fcp 1 , but much less 
than the coherence length of the superconductors. In this limit when kpL 3> 1, one can 
neglect intermode scattering within the constriction and write the electron wave function 
as a superposition of different independent modes for the transverse motion. As a result, 
both current and noise is a sum of independent contributions from each mode. Each term 
corresponding to a certain mode can be expressed through the solution of the scattering 
problem for electronic states in the leads. In the following, a general theoretical approach 
to non-equilibrium noise in voltage-biased contacts is presented and noise is given in terms 
of the solution of the scattering problem. 

The case of perfect transparency, D = 1, is particularly interesting because then there is 
no excess noise in the normal state (see below). In this limit we give an explicit result for 
large bias voltages, eV 3> A. In the superconducting state there is a finite level of excess 
noise due to Andreev reflection. The physical origin is that the amplitude for this type of 
reflection is non-zero but less than one, for quasiparticle states outside the gap, i.e. having 
energies \E\ > A. The quasiparticle then have a finite probability to be reflected from, or 
transmitted through, the junction which is required in order to have a fluctuating current. 
We show that the zero-frequency limit of the noise spectrum is proportional to the excess 
current [|]||, the proportionality coefficient being e/5. 

This letter is organized as follows. First, a formulation of the general scattering problem 
will be given and the solution for D = 1 specified. Then the excess noise will be expressed 
in terms of the scattering solution, and finally we apply this to the limiting case D = 1, zero 
temperature and zero frequency in the large voltage limit. 
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Scattering problem. In the original Landauer picture |JL5|| , scattering in a normal voltage- 
biased point contact is elastic. The voltage V only enters as a relative displacement of the 
Fermi levels in the reservoirs. For a contact between superconductors, an alternative picture 
where the voltage V is included in the phase of the wave function is more convenient. 
According to this picture, constant electric potentials in the leads far from the contact are 
removed by local gauge transformations which are different on the two sides of the junction. 
Consequently, the voltage is included as phase factors in all scattering states. Transmission 
through the contact now formally appears as inelastic and for a normal contact there is one 
transmitted state at energy E + eV with amplitude d, and one reflected state with amplitude 
r at the same energy E as the incoming state. In a voltage-biased superconducting contact 
though, multiple Andreev scattering processes to all energies E + neV (where n is an integer) 
are possible [|TJ. 

In the leads, the electronic states can be described as solutions to the Bogoliubov- 
de Gennes (BdG) equation [H| and have the form of two-component spinors \I/ which 



can be specified by the side of the junction (L or R) as a subscript, and by the di- 
rection of the longitudinal momentum (±) as a superscript. Introducing slowly vary- 
ing spinors the total wave function for each transverse mode will be expanded as 
\& = e - lEt / h J2±e ±tkFX & ± . The spinors $ satisfy a reduced time-independent BdG equa- 
tion ||16|| , (±hkpp/m e + <t x A)$ ± = E$>^ where m e is the electronic mass and p the mo- 
mentum operator. The Fermi wave vector (energy) is kp (Ep) outside the channel, and 



kp = h 1 ^2m e Ep[l — (mr /kpw) 2 } for mode n inside the contact of width w. The spinors $ 

change on the scale of the coherence length £ = h 2 kp/7iAm e , i.e. they can be considered as 
being constant in the contact region. In the following we adopt the Andreev approximation 



|17| which is valid up to lowest order in the small parameter A/E-p. In this approximation 



the scattering is described by the boundary condition |5|.|] . 

( M = v V=( f de-»" vt /* \ ( , 

\®r) \®r) " V de*' eVt/h -r*d/d* )' 1 1 

The parameter r(d) is the normal reflection (transmission) amplitude, while the matrix 
f(d) is the same quantity multiplied by the 2x2 unit matrix. Since the junction length 
is much less than the coherence length |]I8|| , the coordinate dependence of the phase of 



the order parameter can be approximated by a step function centered at the junction, and 
the magnitude of the order parameter by a constant. This allows us to include the phase 
difference eVt/H in the matching matrix V, see above. The general boundary condition (|]) 
has to be applied to the four distinct scattering cases - electron- and hole-like quasi particles 
incident from the left and the right, and for positive as well as for negative energies. 

The scattering problem substantially simplifies for the case of perfect transmission, d = 
1. The wave vector is then conserved to lowest order in A/Ep. Hence only wave vector 
conserving scattering events can take place, that is, Andreev reflection to energies E + 2neV 
with amplitudes Ai n and normal transmission to energies E + (2n + l)eV with amplitudes 
A/2n+i- Further, the direction of energy transfer is prescribed by the sign of the voltage 
and by the scattering case. For a given scattering case only scattering to higher (or lower) 
energies is possible. A typical scattering state is shown in Fig. |]. The scattering amplitudes 
and A/2n+i can be expressed as 
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•A.2n — A 2n , Nln+l ~ A 2n+ i ? (2) 

through the function X n , 

x f [(«o - vl)/v Q u n \ UkZl(v k /u k ), n^O 
n \ -vq/uq, n = . 1 ; 

with upward (downward) recursion, n > 0, for scattering states associated with incident 
electron-like (hole-like) quasiparticles from the left and hole-like (electron-like) ones from 
the right. The coherence factors u n = u(E + neV), v n = v(E + neV) are for E > given by 



w 2 (£) = min[l,£/A](l + ^/l - A 2 /E 2 )/2 and t> 2 (£) = min [1, E/A}(1 - yjl - A 2 /E 2 )/2; 
for E < 0, w(-E') = —v*(—E) and f(-E) = u*(—E). In the large voltage limit eU ^> A only 
three scattering amplitudes have to be taken into account, 

A = -v /u , A± 2 = -v±i/u ±1 , A4i = (wq - Vq)/(u u ±1 ) . (4) 

All other amplitudes are at most of order A/eV for all energies. 

Excess Noise. Current noise is usually defined as the Fourier transform of the current- 
current correlation function, S(r) = l/2( I(t + r),I(t) ) - (I(t)) 2 , where [A, B] + = AB + 

BA. Using the Bogoliubov-de Gennes transformation JjiJ, we can express electron field 
operators in terms of quasiparticle operators 7. The current operator will be 

Hi) = [l\l23+{\%t) + 7I7 ±7_ (1,2,0 + 7172^(1, 2, *) +7i t 7 2 Wi(2,l,t)] • (5) 

The sum is over quantum labels 1 and 2 ranging over all states (meaning scattering case 
index, energy of incoming state, side bands n, and spins a). Above we have defined three 
different types of generalized currents, 

j+(l,2,t) = ($(l)|iV<7n$(2)) + (iV$(I)|<7u$(2)>, (6) 
j_(l,2,t) = ($*(!) |zVa 22 <I>*(2)) + (zW(l)|<7 22 <I>*(2)>, (7) 
4,(1,2,*) = ($*(!) |zV<7 21 <I>(2)) + (zW(I)|«7 21 $(2)> , (8) 

expressed through the 2x2 matrices having element (i, j) equal 1 and zero otherwise. The 
scalar product between two spinors is here defined as (assume $ T = (u, v)), ($(1)|$(2)) = 
u*(l)u(2)+v*(l)v (2). Introducing the normalized BCS density of states per spin for incoming 
quasiparticles, g{E) = [\E\/y / E 2 — A 2 ]Q(\E\ — A), our general expression for noise can be 
written as [|19H 



p2 /-oo 

S(cu,V,T) = -—J2 f(E)[l-f(E + neV±hu;)} 

o7T II j s _ n J-fi 



x g(E)g{E + neV±huj)Trtt {n) {E,E + neV±huj)dE, (9) 
fi (n) (E, E') = Q {n) (E, E')P(E', E) + ]T Q<?\E,E')P U (E' \E) , (10) 



u=l,2 



Q^{E,E ! )^Y,^n +m {E)s\^{a uv )^ m {E'))^ (11) 

m 

P (u) (E', E) = (£ &n(E')s\ ®(<r w ) E * m (E)) . (12) 
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Above, $> m (E) is a row vector for harmonic side-band m consisting of four two-component 
(column) spinors <& m (E), each corresponding to one of the four scattering cases. The symbol 
denotes outer product in scattering case indices. Thus Q and P are 4x4 matrices 
with element corresponding to scattering cases % and j in the bra and ket states in 

the definitions above. The trace is then calculated in this four-dimensional space. The 
parenthesis around a vv implies that <j vv should be excluded when calculating and P. 
The operator s causes the inner product to be multiplied by the average of the wave vectors 
(in units of the Fermi vector) of the bra and ket states - only if the wave vectors have equal 
sign will the term contribute (i =>■ ±1). 

For perfect transmission (d = 1), the wave vector is conserved when quasiparticles are 
scattered in the junction implying that all subbands have the same wave vector (to leading 
order in A/Ep). This allows us to assign a value of the sign of the wave vector (±kp) to 
each scattering state containing all its subband components. The number of contributing 
terms in the trace in the noise expression is then reduced from 16 to 8 (see operator s in 
S(u>)). In the limit of large bias voltages V all scattering amplitudes, to lowest order in 
A/eV, are given by Eq. (f|). 

Putting in proper coherence factors into the expressions for scattering amplitudes and 
spinors, we find that for zero temperature, the zero frequency limit of the noise can be 
written as 

Se X c(0) - H D A (E)[1 - D A {E)\ dE = ^ , (13) 

as A/eV — > 0. The quantity D A (E) = 1 — Aq can be interpreted as the elastic Andreev 
transparency at energy E. As a function of E, D A = for < A (complete Andreev 
reflection). It is close to 1 for E ^> A, while in the vicinity of the gap for \E\ > A, both D A 
and 1 — D A are finite. For arbitrary voltages, noise also contains inelastic contributions. In 
the limit of large voltages though, these terms vanish because of the decay of the Andreev 
reflection probability Aq far away from the gap region \E\ < A. This voltage independent 
level of noise is proportional to the excess current 5exc(0) = |e/ exc |/5. It is illuminating 
to compare this with the level of noise in a normal quantum point contact JTT|-|T3"[ at u = 
and T = 0, 

Sn(0) = ^-[° D N (E)(l-D N (E))dE, (14) 

7T ll J —eV 

where D^{E) is the normal transparency. In our case, Dn = \d\ 2 = 1, and Sn vanishes. 

Conclusion. We have shown that there is a finite level of current noise in voltage-biased 
junctions between superconductors, even if the transparency in the normal state is perfect. 
This excess noise is induced by fractional Andreev reflection and for large bias voltages 
(eV 3> A) it is proportional to the excess current, the proportionality coefficient being equal 
to \e\/5. Consequently, the relationship between the excess noise and the excess current is 
qualitatively similar to the classical relation between shot noise and average current in the 
low- voltage regime, S = 2\e\I. The quantum reduction of excess noise (coefficient 1/5) is in 
our case pronounced. 
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work on the Physics of Nanometer Electronic Devices, the Swedish Natural Science Research 
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FIG. 1. A scattering state for the perfect transmission (d = 1) case. Full lines represents 
one quasiparticle character e(h) while dashed lines means the opposite character h(e) - full thick 
line corresponds to incoming quasiparticle. Some Andreev reflection A n and normal transmission 
amplitudes M n to sidebands n with energy E+neV are shown. Interchanging e and h characters, or 
changing the sign of the applied bias voltage, or interchanging left and right, causes the scattering 
state to be mirror-imaged around the energy E of the incoming particle, n — > — n. 
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